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Abstract 

Background: Pseudomnas putida is a natural producer of medium chain length polyhydroxyalkanoates (mcl-PHA), 
a polymeric precursor of bioplastics. A two-fold increase of mcl-PHA production via inactivation of the glucose 
dehydrogenase gene gcd, limiting the metabolic flux towards side products like gluconate was achieved before. 
Here, we investigated the overproduction of enzymes catalyzing limiting steps of mcl-PHA precursor formation. 

Results: A genome-based in silico model for P. putida KT2440 metabolism was employed to identify potential 
genetic targets to be engineered for the improvement of mcl-PHA production using glucose as sole carbon source. 
Here, overproduction of pyruvate dehydrogenase subunit AcoA in the P. putida KT2440 wild type and the Aged 
mutant strains led to an increase of PHA production. In controlled bioreactor batch fermentations PHA production 
was increased by 33% in the acoA overexpressing wild type and 121% in the acoA overexpressing Aged strain in 
comparison to P. putida KT2440. Overexpression of /^/-encoding 6-phosphoglucolactonase did not influence PHA 
production. Transcriptome analyses of engineered PHA producing P. putida in comparison to its parental strains 
revealed the induction of genes encoding glucose 6-phosphate dehydrogenase and pyruvate dehydrogenase. 
In addition, NADPH seems to be quantitatively consumed for efficient PHA synthesis, since a direct relationship 
between low levels of NADPH and high concentrations of the biopolymer were observed. In contrast, intracellular 
levels of NADH were found increased in PHA producing organisms. 

Conclusion: Production of mcl-PHAs was enhanced in P. putida when grown on glucose via overproduction of a 
pyruvate dehydrogenase subunit (AcoA) in combination with a deletion of the glucose dehydrogenase {gcd) gene 
as predicted by in silico elementary flux mode analysis. 

Keywords: Medium chain length polyhydroxyalkanoate, Pseudomonas putida, Systems metabolic engineering, 
Pyruvate dehydrogenase, Glucose 



Background 

Metabolic engineering of microorganisms has significantly 
contributed to the biotechnological production of chemicals 
in a sustainable manner [1,2], Insertion of heterologous 
pathways and modification of the host metabolism, 
via deletion or amplification of genes, has been a common 
procedure to enhance the synthesis of a desired compound 
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[3,4], Nevertheless, these modifications often evoke 
detrimental effects on the cell metabolism which lead to 
decreased growth and low performance of the constructed 
cell factory [2,3,5]. In silico reconstruction of the metabolism 
of bacterial production hosts in combination with metabolic 
flux modeling can help to design the optimal pathways 
towards a product of interest. Multiple efforts have 
been made to engineer microorganisms to obtain 
sustainable polymers which have the potential to replace 
petroleum-based plastics. In this regard, medium- 
chain-length polyhydroxyalkanoates (mcl-PHAs) are among 
the most promising biodegradable polymers synthesized 
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by microorganisms which can be used in a wide range of 
applications [6]. Bacteria from the genus Pseudomonas are 
natural producers of mcl-PHAs since they have the entire 
machinery to synthesize these polyesters from different 
carbon sources (Figure 1) [7]. They accumulate these 
polymers as inclusion bodies in the cytoplasm promoted 
preferentially by a high carbon supply and the limitation 
of an inorganic nutrient such as nitrogen, oxygen, or 
phosphorous [8]. As other PHA-accumulating bacteria, 



Pseudomonas putida strains utilize PHAs as a reservoir 
of carbon and energy to cope with the changing 
environmental conditions of their natural habitats. 
For decades, studies for biotechnological mcl-PHA 
production focused on the use of fatty acid as carbon 
source [9-12]. This approach led to PHA accumulation in 
the cell, yielding more than 80% of its dry weight as PHA 
[13]. Additionally, engineering of the (3-oxidation pathway 
in P. putida allowed for production of novel mcl-PHAs, 
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Figure 1 Central metabolism of P. putida KT2440. Enzymes of interest are indicated. Enzymatic steps with NADVNADP and NADPVNADPH 
involvement are marked. 
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which implies a broad range of applications in biomedicine 
and for biomaterials [10,14,15]. Much less attention has 
been given to the engineering of metabolic pathways 
involved in sugar or polyalcohol utilization in P. putida 
[16]. It was reported that inactivation of the isocitrate lyase 
gene resulted in a 1.7-fold increased production of 
mcl-PHA using gluconate as carbon source [17]. 

Different waste materials were tested as feedstock 
for PHA producing Pseudomonas species in order to 
reduce the overall production cost. Pseudomonas spe- 
cies accumulated mcl-PHAs when grown on sodium 
terephthalate produced from a PET pyrolysis, raw gly- 
cerol [18], polystyrene pyrolysis oil, and animal waste 
lipids [19]. Glucose is still one of the cheapest and most 
available feedstock for industrial polymer biosynthesis 
[20]. However, only small amounts of mcl-PHAs are 
currently synthesized utilizing glucose as carbon source 
[21,22]. Most metabolic approaches for the improvement 
of PHA production with glucose as carbon substrate were 
performed using Escherichia coli [23-25]. After cloning 
the PHA synthesis genes from Cupriavidus necator 
(formerly known as Ralstonia eutropha) into E. coli higher 
levels of NADPH versus NADP + were found essential to 
promote PHB synthesis. Consequently, bioengineering 
efforts on the rising of NADPH availability in turn to 
improve PHA synthesis on glucose were undertaken. 
Overproduction of glucose-6-phosphate dehydrogenase 
(G6PDH encoded by zwf) increased the intracellular 
NADPH levels and led to an increased accumulation of 
PHA by 41% [26]. A similar approach was carried out with 
C. necator, which resulted in an increased production of 
PHA in the presence of fructose [27]. It is known that 
Pseudomonas species prefer other carbon sources like 
organic acids and amino acids over glucose [28]. 
Consequently, our research focused on the development 
of P. putida strains which can accumulate high amounts 
of mcl-PHA on glucose as the sole carbon source. 
Recently, novel bioengineering approaches for an 
optimized production from glucose were undertaken 
[21]. Initially, a genome-based in silico model of the 
P. putida metabolism was established and used for 
the prediction of optimal fluxes towards PHA and its 
precursor molecules. Among a list of different genetic 
targets, elementary flux mode analysis suggested an 
inactivation of glucose dehydrogenase encoded by gcd to 
prevent undesired by-product formation and excretion. In 
agreement with the prediction from the in silico model 
the constructed P. putida Aged mutant strain showed an 
increased PHA production in batch [13] and Fed-Batch 
cultivations [29]. Here, we continued to experimen- 
tally verify the predictions of our in silico model for 
the optimization of mcl-PHA production in P. putida. 
We tested for the proposed positive impact of the 
overexpression of the genes for pyruvate dehydrogenase 



(acoA) and 6-phosphogluclactonase (pgl) on mcl-PHA 
production. 

Results and discussion 

Overexpression of the genes for the pyruvate 
dehydrogenase AcoA and 6-phosphoglucolactonase Pgl 
for improved PHA production 

The acoA and pgl genes were separately cloned into the 
broad host range vector pSEVA424 and then transformed 
into the P. putida KT2440 and the Aged mutant 
strains. A list of all used strains, vectors and primers 
are summarized in Table 1. Wild type P. putida KT2440 
and the corresponding Aged mutant strain carrying the 
vector without insert served as controls. Cells were grown 
in defined M9 mineral minimal medium supplemented 
with 20 (g/L) glucose and 62.5 [ig/mL streptomycin as 
described in the Methods section. Protein production 
was induced at an initial OD 600nm of 0.05. Cells were 
harvested after 60 hours by centrifugation prior to 
protein composition analyses and PHA-extraction and 
quantification. SDS-PAGE analysis was used to assess 
the overproduction of the AcoA and the Pgl proteins. 
Successful overproduction of AcoA (M r = 35,000) and Pgl 
(M r = 23,000) was observed, which is in good agreement 
to the deduced relative molecular mass of approximately 
34,700 Da for AcoA and 25,500 Da for Pgl (Figure 2). 

Increased mcl-PHA production in P. putida by overexpression 
of the pyruvate dehydrogenase gene acoA 

In order to evaluate the effect of acoA and pgl overex- 
pression on PHA versus biomass production GC-MS 
analyses of the extracted PHA and biomass determination 
were performed (Figure 3). Triplicate aerobic cultures of 
P. putida strains KT2440 (wt), KT2440-^coA, KT2440-^/, 
KT2440A^, KT2440Agcd-acoA, and KT2440 Agcd-pgl 
were tested. For this purpose, cells were grown in defined 
minimal medium for 60 hours until glucose as the 
sole carbon source was completely consumed (in all 
tested strains no glucose was detected via HPLC at 
60 hours). As predicted by the in silico design the 
pyruvate dehydrogenase-producing KT2440-^coA and 
KT2^0 Aged- aco A mutant strains both showed an in- 
creased PHA content amassing 33.3 wt % and 42.1 wt 
% as mcl-PHA, respectively (Figure 3A). This was not 
the case for the phosphoglucolactonase-overproducing 
KT2440-^/ (25.7 wt %) and KT2440Agcd-pgl (33.2 wt %) 
strains in which the PHA concentration remained at the 
basic level. In addition, no obvious major changes of the 
mcl-PHA momomer composition between the various 
producing P. putida strains were observed (Table 2). Next, 
the yield of biomass of engineered strains grown on 
glucose was determined. As initially desired, the biomass 
production of the various recombinant strains was 
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Table 1 Bacterial strains, plasmids and primers used in this study 



Strains 



Relevant features/genotype* 



Source reference 



E. coli 

DH10b 

XL10-Gold 

DH10b-6?co/\ 

DH10b-pg/ 

XL10-Gold- acoA 

XL10-Gold-pg/ 

P. putida 

KT2440 

KY2440-ocoA 

KX2440-pg/ 

KX2440Agcc/ 

KY2440Agcd- 
ocoA 

KY2440Agcd-pgl 
Plasmids 

pJETI .2 

pSEVA424 

pM].2-acoA 

pMM-pgl 

pSEVA- acoA 

pSEVA-pg/ 

Primers 

acoA Fw 

acoAWv 

pg/Fw 

pg/Rv 



F" mcrA Mmrr-hsdRMS-mcrBQ O80d/acZAM15 AlacX74endA]recA]deoRA{ara,leu) 
7697araD]39galVgalKnupGrpsL X" 

end A] glrN44recA} thi-1 gyr A96relAl loc\-\XeA{mcrA) 1 83 A(mcrCB-hsdSMR-mrr) 1 73tet R 
F'[proAB/ac/ q ZAM15 Tn70(Tet R Amy Cm R )] 

DH10b carrying the P. putido ocoA gene cloned in pJETI .2 

DH10b carrying intermediate pgl gene cloned in pJETI .2 

XL10-Gold harboring acoA overexpression vector pSEVAaco 

XL10-Gold harboring pgl overexpression vector pSEVApg/ 

Wild type 

KT2440 harboring ocoA expression vector pSEVA424 
1^2440 harboring pgl expression vector pSEVA424 
1^2440 carrying gcd mutation 

1^2440 carrying the gcd knockout mutation harboring acoA overexpression vector pSEVAc/co/\ 
KT2440 carrying the gcd knockout mutation harboring pgl overexpression vector pSEVApg/ 

Carrier vector, toxin expressed upon self-ligation, Ap r 
Expression vector for P. putido, IPTG-inducible, toe promoter, Sm r 
Plasmid pJETI .2 carrying ocoA gene for propagation, Ap r 
Plasmid pJETI .2 carrying pgl gene for propagation, Ap r 
ocoA overexpression vector, based on pSEVA424, Sm r 
pgl overexpression vector, based on pSEVA424, Sm r 

5 CCGGAATTCJKGGJKGGAAAAACA 7ATGTCCA ATC A ACTCAGT3 

5TGC7~C77\G/\GCATCAGGGGTAGGCGACGTA3' 

5'CCGGA47TCAGGAGGAAAAACflrATGGGAGGGCGTGGTATG3 

5TGC7C77\G/\GCATCATGGGCACCAGTAGAT3' 



Invitrogen 

Stratagene 

This study 
This study 
This study 
This study 

DSMZ 

This study 

This study 

Poblete et al., 2013 
[21] 

This study 
This study 

Fermentas 
De Lorenzo's Lab 
This study 
This study 
This study 
This study 

This study 
This study 
This study 
This study 



Ap r = ampicillin resistance; Sm r = streptomycin resistance. 

¥ Source the material was obtained from or reference. EcoR\ and Xba\ restriction sites are shown in italics. Shine Dalgarno (RBS) sequences are displayed in bold 
whereas the eight conserved nucleotides upstream the RBS are depicted in bold and italics. 



found very similar to the one determined for the wild 
type P. putida KT2440 strain (Figure 3B). 

In summary, the most efficient PHA producer was 
KT2440Agcd-acoA strain, where the PHA concentration 
was doubled compared to P. putida wild type strain 
(Figure 3 A) and enhanced by 21% compared to the 
KT2440Agcrf mutant strain (Figure 3A). Thus, the 
predicted engineering strategy was successful in in- 
creasing the mcl-PHA synthesis of P. putida. 

Batch fermentation in the bioreactors for controlled 
mcl-PHA production 

To ensure well-controlled mcl-PHA production using 
our metabolically engineered P. putida strains, aer- 
obic batch fermentations of the KT2440-acoA and 



KT24 ! 4 ! 0Agcd-acoA strains on glucose were performed in a 
bioreactor. Figure 4 shows the time-resolved measure- 
ments for mcl-PHA formation, biomass production and 
the corresponding concentrations of ammonium, gluco- 
nate, 2-ketogluconate, and glucose formation. As expected, 
the KT2440-acoA strain, gluconate and 2-ketogluconate 
were excreted into the medium during consumption of the 
available glucose (Figure 4 A, Table 3). In contrast, during 
the fermentation of the flux optimized KT2440 Aged- acoA 
strain (Figure 4B) excreted organic acids were below the de- 
tection limit (Table 3). A recent study compared mcl-PHA 
synthesis in P. putida KT2440 and KT2442 when grown on 
gluconate [30], P. putida KT2440 accumulated 17% of its 
CDW as PHA, whereas KT2442 only synthesized 1.7 wt % 
as PHA, demonstrating the large differences between the 
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Figure 2 SDS-PAGE gels demonstrating the recombinant production of pyruvate dehydrogenase and 6-phosphogluconolactonase 
by the recombinant P. putida strains. (A) Shown the separations of whole cell protein extracts from 1^2440 (lane 1), \<J2440-ocoA (lane 2), 
r^2440-pg/ (lane 3), and (B) KY2440Agcd (lane 4), r^2440Agcd-aco/\ (lane 5) and r^2440Agcd-pg/ (lane 6). 



two P. putida strains. P. putida KT2440 can re-consume 
gluconate and 2-ketogluconate once secreted to the 
medium. We recently demonstrated that this process is not 
efficient to improve mcl-PHA synthesis [29], resulting in a 
lower volumetric PHA productivity as compared to the one 
shown by the Aged P. putida mutant strain under the same 
growth conditions. Several investigations have shown 
that the volumetric productivity is the key parameter for 
the industrial production of PHAs and thus, any factor 
influencing this parameter affects the economics and 



cost competitiveness of PHAs synthesis via microbial 
fermentation [30-33]. Ammonium was completely con- 
sumed at 12 hours by the KT2440-<zcoA strain, while 
KT2440Agcd-acoA took more than 15 hours to deplete 
the entire amount of nitrogen. In addition, a large propor- 
tion of carbon was directed towards gluconate in the 
KT2440-<2coA strain. More than 5 g/L of gluconate was 
produced after 24 h, (Figure 4A). Furthermore, the lower 
biomass of 2 g/L obtained for this strain prior ammonium 
limitation was most likely due to the increased production 
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Figure 3 mcl-PHA (A) content and cell dry weight (B) of different metabolically engineered indicated P. putida strains. 



of cellular compounds including nucleic acids and 
proteins due to less PHA formation. These findings might 
explain the reduced PHA and biomass yield on glucose in 
comparison to the KT 2440 Aged- a coA mutant strain. 
The genetic modification of the KT244!0Agcd-acoA led to 
the lowest u max of all strains (0.25 h" 1 versus 0.56 h" 1 for 
KT2440). Nevertheless, the KT2440 Aged- acoA strain 
biomass (CDW = 4.71 g/L after 50 hrs) and mcl-PHA 
42.1 wt % contents were higher in comparison to the 
KT2AAQ-acoA strain (CDW = 3.63 g/L after 50 hrs and 
mcl-PHA = 33.3 wt %). To our knowledge this is the high- 
est mcl-PHA accumulation ever reported using a P. putida 
strain in the presence of glucose in batch cultures. 

Transcriptome analysis of the engineered P. putida PHA 
production strains 

Glucose is metabolized in P. putida that utilizes three path- 
ways converging at the level of 6-phosphogluconate prior 
to the formation of 2-keto-3-deoxy-6-phosphogluconate 
(KDPG; [34]). The repressor HexR regulates the expression 
of major genes of the pathway including the zwf-1, pgl, eda 
operon and the gapl gene encoding glyceraldehyde-3- 
phosphate dehydrogenase. Accumulation and binding of 



KDPG to HexR lead to derepression [35-37]. Furthermore, 
the genes of PHA synthesis are subjected to a complex 
control by catabolite regulation and various stress re- 
sponses mediated by the complex RNA-dependent GacS 
and Crc regulating systems, besides others [38-40]. In 
order to understand the gene regulatory adaptation of the 
newly engineered strains to enhanced PHA produc- 
tion the gene expression compared to the parental 
strains was analyzed using next generation sequencing 
of transcripts (RNA Seq ). We are fully aware that major 
adaptation processes might take place at the level of 
enzyme regulation which cannot be made visible with 
the employed approach. Samples were taken after 
60 hours of cultivation, when the PHA concentration in 
the cell reached its highest level. Strict cut-off parameters 
(log change > 2 and a j^-value < 0.05) were applied. First, 
the gene expression of acoA overexpressing strains 
was compared to the wild type KT2440 strain and the 
KT2440Agrd mutant strain (Table 2). Overall, 53 open 
reading frames (ORFs) were found differentially tran- 
scribed (12 up-regulated and 41 down-regulated) between 
the KT2440 wild type strain with or without overex- 
pressed acoA. Moreover, 108 ORFs (59 up-regulated and 



Table 2 Monomer composition of mcl-PHA produced by metabolically engineered P. putida strains in flask experiments 



Strain C6°° C8°° C10 00 C12°° CI 2:1 00 C14 c 



KT2440 


n.d. 


1 1 .8 ± 0.2 


73.0 ±0.1 


4.7 ± 0.2 


9.4 ± 0.7 


0.8 ± 0.3 


KT2440Agcc/ 


n.d. 


14.0 ±0.4 


71.1 ±0.8 


4.6 ± 0.9 


8.5 ± 0.4 


1.1 ±0.2 


KT2440-acoA 


n.d. 


13.0 ±0.6 


69.1 ±1.2 


3.4 ± 0.7 


10.5 ±0.2 


3.9 ± 0.2 


KT2440bgcd-acoA 


n.d. 


14.8 ±0.7 


70.2 ± 1 .4 


5.1 ±0.3 


8.8 ±0.1 


0.7 ± 0.3 


KT2440-pgl 


n.d. 


12.2 ±0.9 


74.1 ±1.5 


3.8 ± 0.5 


9.9 ± 0.8 


n.d. 


KT2440Agcc/-pg/ 


n.d. 


14.6 ±0.3 


73.1 ±0.8 


3.7 ± 0.2 


8.1 ±0.5 


0.5 ±0.1 



n.d.: not detected, less than 0.2%. 

°°: The data were determined by GC/MS and are given as relative molar fraction (%) of C6: 3-hydroxyexanoate, C8: 3-hydroxyoctanoate, C10:3-hydroxydecanoate, 
C12: 3-hydroxydodecanoate, CI 2:1: 3-hydroxy-5-cis-dodecanoate, and C14: 3-hydroxytetradecanoate. 
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Figure 4 Growth and mcl-PHA production of the metabolically engineered strains KY2440-acoA (A) and KT2440Ago/-acoA (B) in batch 
cultures. Cells were grown in lab scale bioreactor in the presence of glucose 18.5 (g/L). Data are presented as mean values obtained from two 
independent experiments. 



49 down-regulated) were found changed in their transcrip- 
tome between the KT2440Agr<i strain with or without 
overexpressed acoA. As expected, the strongest up- 
regulated gene in both variant strains overexpressing 
acoA was PP_0555 (acoA). Interestingly, also PP_0554 



(acoB) was found upregulated, however to a lesser extent 
(Table 4). Most of the genes encoding for enzymes 
related to the central and PHA synthesis metabolic 
pathways showed only slight variations in their expression 
levels (Table 4). Interestingly, the increased flux toward 
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Table 3 Physiological parameters of the wild type P. putida KT2440 and metabolically engineered strains in batch 
culture grown with of glucose 



Strain 


Mmax (h ^ 


Yx/s § (g/g) 


Y2-ketoglut/S (g/g) 


Yg|uconate/S (g/g) 


y pha/s (g/g) 


KT2440 


0.56 


0.22 


<0.001 


0.64 


0.05 


KT2440Agccf 


0.42 


0.23 


<0.001 


<0.001 


0.08 


KY2440-acoA 


0.40 


0.21 


0.04 


0.38 


0.07 


KT2440bgcd-acoA 


0.25 


0.25 


<0.001 


<0.001 


0.11 


KT2440-pg/ 


0.33 


0.21 


0.07 


0.55 


0.05 


KT2440Agcc/-pg/ 


0.21 


0.24 


<0.001 


<0.001 


0.08 



§ Y x/S is calculated based on the total biomass concentration (true biomass + PHA content). Standard deviations of less than 10% were obtained in 3 
independent experiments. 



acetyl-CoA and as a consequence citrate induced the 
expression of the isocitrate lyase gene, which encodes 
the first enzyme of the glyoxylate shunt. One metabolic 
consequence is a reduced NADH production in the TCA 
cycle. A second consequence might be the channeling of 
acetyl-CoA into the glyoxylate cycle. In agreement, it was 
previously shown that inactivation of the isocitrate lyase 
gene leads to increased PHA production [17]. Overall, 
the gene expression data do not provide conclusive 
explanations for the observed increase in PHA produc- 
tion. Most likely regulatory effects at the enzyme level 
might be responsible. 

Next, the pgl overexpressing strains without increased 
PHA production were analyzed for potential compensatory 
gene regulatory activities. Of the 86 ORFs differentially 
expressed in KT2440-pgl mutant strain compared to the 
wild type strain, 77 corresponded to genes that were tran- 
scriptionally up-regulated, whereas 9 genes were found with 
decreased expression levels (Table 5). For KT2440Agcd-pgl 
mutant strain compared to the KT2440Agcd mutant, 15 
genes were found down-regulated while 117 ORFs were 
found induced. As expected, the pgl overexpressing strains 
revealed an increase of pgl mRNA. Glucose 6-phosphate 
dehydrogenase (encoded by PP_1022, zwf) was found to be 
highly up-regulated (>12 fold) in both strains harboring the 
overexpression system (Table 5). Obviously, the initial goal 
to enhance the flux from glucose toward KDPG was 
achieved. The increased amounts of the metabolite might 
have enhanced zwf-1 expression via HexR derepression. 
Furthermore, expression differences were detected for the 
pyruvate metabolism, where the genes encoding pyruvate 
dehydrogenase (PP_0554), pyruvate carboxylase (PP_5347) 
and phosphoenoylpyruvate (PP_1505) showed differential 
expression levels (Table 5). At least for the wild type 
strain harboring the overexpressed pgl a problematic 
rearrangement of pyruvate metabolization occurred 
which might be responsible for the observed production 
phenotype. While the pyruvate dehydrogenase gene acoA 
was found reduced, the pyruvate carboxylase gene acoC-2 
was found induced (Table 5). In addition, induction of the 
phosphoenolpyruvate carboxylase gene ppc was detected. 



The KT2440A < gr<i mutant carrying overexpressed pgl did 
not reveal this response. For the pgl overexpression 
the observed gene regulatory response provides reasonable 
explanations for the observed PHA production phenotype. 

The role of NADPH and NADH in P. putida PHA 
production 

The production of mcl-PHAs in P. putida KT2440 
requires NADPH + H + . P. putida catabolizes glucose via the 
Entner-Doudoroff pathway and hence the NADPH + H + 
producing glucose-6-phosphate dehydrogenase (Zwf) 
is highly active. However, the attempts to further increase 
the metabolic flux through this step by overproducing the 
rate-limiting, succeeding enzyme (Pgl) failed to increase 
PHA production (Figure 1). However, the flux of carbon 
towards PHA was most likely hampered at the pyruvate 
metabolization step in the pgl overexpression strain. 
Interestingly, the transcriptome analyses revealed an 
increased zwf gene expression. In order to determine 
the contribution of the NADPH + H + level in P. putida to 
PHA production, NADPH + H + versus NADP + levels were 
determined for the various bioengineered strains (Table 6). 
The ratios of NADPH to NADP + were measured using 
UV/Vis spectrometry and quantified according to their 
colorimetric changes at A = 450 nm. At the maximum 
mcl-PHA accumulation in the cell (60 h culture), we found 
the NADPH/NADP + ratios of all the PHA-overproducing 
strains were lower compared to the corresponding wild 
type P. putida KT2440 (Table 6). Clearly, the employed 
metabolic engineering strategy increased the polyester 
production in both strains (Table 3, Figure 3). As a 
consequence of this increased flux of carbon into PHA 
synthesis, more NADPH is needed to turn 3-ketoacyl-ACP 
into (R)-3-hydroxyacyl-ACP (Figure 1). Therefore, the 
levels of NADPH/NADP + were lower for the engineered P. 
putida strains (Table 6), in which pyruvate dehydrogenase 
was overproduced. A totally different picture was observed 
for the strains in which the pgl gene was overexpressed. 
Here, an increase of the NADPH/NADP + levels was de- 
tected (Table 6). As already deduced from the transcrip- 
tome data, the initial flux through the Entner-Doudoroff 
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Table 4 Expression profile of the genes belonging to PHA biosynthesis and central metabolic pathways in the 
metabolically engineered P. putida strains compared to their parental strains 



Gene name 



Locus tag 



Description 



Fold change 
Agcd-acoA 



KT2440-acoA 



PHA synthesis 

phol PP5008 

phoF PP5007 

phoO PP5003 

phoC2 PP5005 

phaZ PP5004 

phaD PP5006 

phoG PP1408 
Transporters 

oprB-1 PP1019 

oprB-2 PP1445 
Glycolysis/gluconeogenesis 

glk PP1011 

pgi PP1808 

fbp PP5040 

fda PP4960 

tpiA PP4715 

gap) PP1009 

gap2 PP2149 

pg/c PP4963 

pgm PP5056 

eno PP1612 

pyk PP1362 
Pentose phosphate pathways 

zwfl PP1022 

zw/2 PP4042 

zw/3 PP5351 

pgl PP1023 

gnd PP4043 

gnuK PP3416 

^L//C PP3378 

/cguD PP3376 

rpiA PP5150 

rpe PP0415 

tktA PP4965 

to/ PP2168 
Entner-Doudoroff pathway 

edd PP1010 

eda PP1024 
Pyruvate metabolism 

ocoA PP0555 

acoS PP0554 



PHA granule-associated 
PHA granule-associated 
PHA polymerase 
PHA polymerase 
PHA depolymerase 
Transcriptional regulator 
Acyl-transferase 

Porin 
Porin 

Glucokinase 

Glucose-6-phosphate isomerase 
Fructose-1 ,6-bisphosphatase 
Fructose-1 ,6-bisphosphate aldolase 
Triosephosphate isomerase 
GAP dehydrogenase, type I 
GAP dehydrogenase, type II 
Phosphoglycerate kinase 
Phosphoglyceromutase 
Phosphopyruvate hydratase 
Pyruvate kinase 

G6P dehydrogenase 



6-phosphogluconate dehydrogenase 
6-phosphoglucolacto dehydrogenase 
Carbohydrate kinase 
Dehydroglucokinase 
2-Ketogluconate 6-phosphate reductase 
Ribose-5-phosphate isomerase A 
Ribulose-phosphate 3-epimerase 
Transketolase 
Transaldolase B 

6-Phosphogluconate dehydratase 
KDPG aldolase 

Pyruvate dehydrogenase 
Pyruvate dehydrogenase 



0.7 
-0.3 
-0.2 
-0.6 

0.3 
-0.4 

1.5 

-0.7 
-0.5 

0.2 
-0.5 
0.6 
0.9 
0.8 
-1.2 
0.2 
0.2 
0.4 
0.3 
-0.3 

0.9 
0.1 
0.2 
0.5 
-0.6 
-1.1 
0.45 
-0.7 
0.3 
-0.5 
0.3 
0.4 

-0.3 
0.0 

12.4 
1.5 



0.9 
0.2 
0.7 

-0.2 
0.5 
0.3 

-0.1 

-1.0 

-0.6 

-0.2 
0.0 
0.4 
0.1 
0.6 

-0.4 
0.3 
0.1 
0.1 
0.0 
0.4 

0.2 
-0.3 

0.4 
-0.9 
-0.1 
-0.6 
-0.4 
-0.5 

0.0 
-0.1 

0.7 

0.3 

-0.3 
-0.5 

5.3 
1.0 
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Table 4 Expression profile of the genes belonging to PHA biosynthesis and central metabolic pathways in the 
metabolically engineered P. putida strains compared to their parental strains (Continued) 



acoC 


PP0553 


Pyruvate dehydrogenase 


0.7 


0.7 




PP0545 


Aldehyde dehydrogenase 


0.1 


-1.0 


acsA 


PP4487 


Acetyl-CoA synthetase 


0.8 


1.4 


occC-2 


PP5347 


Pyruvate carboxylase 


0.5 


0.2 


ppsA 


PP2082 


Phosphoenolpyruvate synthase 


0.5 


0.4 


ppc 


PP1505 


Phosphoenolpyruvate carboxylase 


0.5 


0.4 


TCA cycle 










gltA 


PP4194 


Citrate synthase 


0.7 


0.7 


acnA 


PP2112 


Aconitate hydratase 


0.4 


-0.4 


ocnB 


PP2339 


Aconitate hydratase 


0.7 


0.3 


led 


PP401 1 


Isocitrate dehydrogenase 


0.3 


-0.5 


sucA 


PP4189 


2-Oxoglutarate dehydrogenase 


-0.7 


-0.2 


sucD 


PP4185 


Succinyl-CoA synthetase sub alpha 


-0.1 


-0.3 


sucC 


PP4186 


Succinyl-CoA synthetase sub beta 


0.2 


0.3 


sdhA 


PP4191 


Succinate dehydrogenase 


0.3 


0.4 


fumC 


PP0944 


Fumarate hydratase 


0.5 


-1.1 


mdh 


PP0654 


Malate dehydrogenase 


-0.4 


-1.3 


Glyoxylate shunt 










aceA 


PP4116 


Isocitrate lyase 


2.9 


1.3 


glcB 


PP0356 


Malate synthase 


0.8 


0.8 



§ Bold numbers represent a differentiated gene expression pattern. 



pathway might have been induced, leading to a 
KDPG accumulation and an increased zwf-1 expression. 
Consequently, sufficient NADPH was produced. However, 
the increased flux of carbon towards PHA was lost at least 
in one case at the pgl overexpression strains. Overall, gene 
and most likely enzyme activity regulatory phenomena 
prohibited further flux towards PHA production and 
resulted in increased NADPH + H + levels. 

By overproducing the AcoA subunit of the pyruvate 
dehydrogenase complex, there was a significant change on 
NADH/NAD + values in comparison to those found in their 
parental strains (Table 6). P. putida strains with high PHA 
production revealed high NADH/NAD + ratio values. It was 
shown before that intracellular NADH concentrations are 
crucial for cell growth [41,42], Similar observations of the 
importance of the NADH/NAD + ratio for PHA synthesis 
from fatty acids in P. putida were published before [43-45]. 

Conclusion 

We successfully continued our in silico metabolic 
modeling based approach for the improvement of PHA 
formation in P. putida. The predicted overexpression of 
the pyruvate dehydrogenase subunit gene acoA in com- 
bination with the deletion of the glucose dehydrogenase 
gene gcd resulted in an increase of PHA production by 
120%. Potential explanation for the observed PHA 



production phenotypes were derived from gene expression 
and cofactor quantification analyses. Beyond this work on 
microbial biopolymer production, in silico based meta- 
bolic engineering has proven truly valuable for succinate 
[46], 1,4 butandiol [47], bio-ethanol [48], and amino acids 
[49,50], suggesting this strategy as highly promising to 
breed superior producers industrial biotechnology. 

Methods 

Bacterial strains 

Pseudomonas putida KT2440 (DSMZ, Braunschweig, 
Germany) as wild type strain and the gene knock-out 
mutants P. putida Aged [21] were used in this study 
(Table 1). E. coli strain DHlOb (Invitrogen, Darmstadt, 
Germany) was employed for the amplification of 
pJET1.2 derivatives (Promega, Darmstadt, Germany). 
E. coli XLIO-Gold (Stratagene, Waldbronn, Germany) was 
used for replication of the pSEVA plasmids (pSEVA424, 
pSEVA-pgl and pSEV A-acoA). The pSEVA plasmids 
constitute a series of low-copy plasmids for the isopropyl 
1-thio-p-D-galactopyranoside (IPTG) inducible overexpres- 
sion of target genes [51]. Transformation of the P. putida 
strains was accomplished via electroporation applying a 
voltage of 2.5 V, 200 fl and 25 \iF, while Ecoli 
DHlOb was transformed as previously described [52]. 
P. putida KT2440 and its Aged mutant were transformed 
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Table 5 Expression profile of the 
metabolically engineered strains 



genes belonging to PHA biosynthesis and central metabolic pathways in the 
compared to their parental strains 



Gene name 



Locus tag 



Description 



Fold change 
Agcd-pgl 



KT2440-pg/ 



PHA synthesis 

phol PP5008 

phoF PP5007 

phoO PP5003 

phoC2 PP5005 

phaZ PP5004 

phaD PP5006 

phoG PP1408 
Transporters 

oprB-1 PP1019 

oprB-2 PP1445 
Glycolysis/gluconeogenesis 

glk PP1011 

pgi PP1808 

fbp PP5040 

fda PP4960 

tpiA PP4715 

gap) PP1009 

gap2 PP2149 

pgk PP4963 

pgm PP5056 

eno PP1612 

py/c PP1362 

Pentose phosphate pathways 

zwfl PP1022 

zwf2 PP4042 

zw/3 PP5351 

pg/ PP1023 

gnd PP4043 

gnu/( PP3416 

kguK PP3378 

/cguD PP3376 

rp/7\ PP5150 

rpe PP0415 

t/ct/\ PP4965 

tal PP2168 
Entner-Doudoroff pathway 

edd PP1010 

edo PP1024 
Pyruvate metabolism 

ocoA PP0555 

acoS PP0554 



PHA granule-associated 
PHA granule-associated 
PHA polymerase 
PHA polymerase 
PHA depolymerase 
Transcriptional regulator 
Acyl-transferase 

Porin 
Porin 

Glucokinase 

Glucose-6-phosphate isomerase 
Fructose-1 ,6-bisphosphatase 
Fructose-1 ,6-bisphosphate aldolase 
Triosephosphate isomerase 
GAP dehydrogenase, type I 
GAP dehydrogenase, type II 
Phosphoglycerate kinase 
Phosphoglyceromutase 
Phosphopyruvate hydratase 
Pyruvate kinase 

G6P dehydrogenase 



6-phosphogluconate dehydrogenase 
6-phosphoglucolacto dehydrogenase 
Carbohydrate kinase 
Dehydroglucokinase 
2-Ketogluconate 6-phosphate reductase 
Ribose-5-phosphate isomerase A 
Ribulose-phosphate 3-epimerase 
Transketolase 
Transaldolase B 

6-Phosphogluconate dehydratase 
KDPG aldolase 

Pyruvate dehydrogenase 
Pyruvate dehydrogenase 



0.6 
0.3 
-0.8 
-0.7 
-0.4 
0.2 
2.1 

0.4 
-0.4 

0.3 
-0.4 
0.3 
0.6 
0.4 
-0.4 
-0.3 
0.1 
0.5 
0.5 
-0.5 

5.2 

0.2 
-0.1 
9.3 

0.2 
-0.1 
0.6 
0.1 
0.3 
0.2 
0.7 
0.2 

0.2 
0.4 

1.3 

0.7 



-0.3 
-0.4 
-0.1 
-0.4 
0.0 
0.1 
0.2 

-0.1 
-0.4 

0.2 
0.0 
-0.1 
-0.1 
0.5 
-0.2 
-0.2 
0.0 
-0.3 
0.1 
0.6 

3.6 

-0.6 
-0.3 
7.8 

-0.3 
0.3 
0.5 
-0.1 
-0.4 
-0.1 
0.3 
0.1 

0.1 
-0.2 

-2.7 

-0.8 
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Table 5 Expression profile of the genes belonging to PHA biosynthesis and central metabolic pathways in the 
metabolically engineered strains compared to their parental strains (Continued) 



acoC 




PP0553 


Pyruvate dehydrogenase 


0.5 


-1.2 






PP0545 


Aldehyde dehydrogenase 


0.5 


-0.9 


acsA 




PP4487 


Acetyl-CoA synthetase 


-0.3 


-0.6 


occC-2 




PP5347 


Pyruvate carboxylase 


0.8 


2.3 


ppsA 




PP2082 


Phosphoenolpyruvate synthase 


0.2 


-0.1 


ppc 




PP1505 


Phosphoenolpyruvate carboxylase 


0.0 


1.4 


TCA cycle 












gltA 




PP4194 


Citrate synthase 


-0.2 


0.7 


acnA 




PP2112 


Aconitate hydratase 


0.8 


-0.5 


ocnB 




PP2339 


Aconitate hydratase 


-0.2 


0.3 


led 




PP4011 


Isocitrate dehydrogenase 


-0.3 


-0.3 


sucA 




PP4189 


2-Oxoglutarate dehydrogenase 


-0.6 


0.2 


sucD 




PP4185 


Succinyl-CoA synthetase sub alpha 


-0.2 


0.5 


sucC 




PP4186 


Succinyl-CoA synthetase sub beta 


-0.1 


0.1 


sdhA 




PP4191 


Succinate dehydrogenase 


-0.7 


0.3 


fumC 




PP0944 


Fumarate hydratase 


1.1 


-0.4 


mdh 




PP0654 


Malate dehydrogenase 


0.8 


-0.8 


Glyoxylate shunt 












aceA 




PP4116 


Isocitrate lyase 


1.2 


-1.1 


glcB 




PP0356 


Malate synthase 


-0.4 


-0.1 


p-value < 0.05. 

§ Bold numbers represent a differentiated gene expression pattern. 








with the plasmids pSEVA-acoA and pSEVA-pgl encoding primers (acoA/Rv, pgl/Rv) encoded an Xbal restriction site 
the pyruvate dehydrogenase subunit A (AcoA) and and the last 18 nucleotides of target gene sequences 
6-phosphogluconolactonase (Pgl), respectively. (Table 1). Amplified genes were cloned into pJET1.2 

resulting in p]ET1.2-acoA and p]ET1.2-pgl, respectively. 
Genetic techniques and construction of bacterial strains The cloned genes were excised from their host vectors via 
The 978 bp acoA and 729 bp pgl genes were PCR amplified EcoRl and Xbal restriction and subsequently cloned into 
using primer pairs acoAFw/acoARv and pglFw/pglRv, the EcoRl-Xbal sites of pSEVA424 resulting in pSEVA-pgl 
respectively. Used primers sequences are listed in Table 1. and pSEV A-acoA, respectively. 
The forward primers (acoA/Ew, pgl/Vw) of both genes 

encoded an EcoRl restriction site, a consensus Ribosomal Shake flask cultivation 

Binding Site (RBS) sequence (AGGAGG), eight con- P. putida strains were kept as frozen stock in 25% 
served nucleotides downstream (AAAAACAT) and 18 glycerol at -80°C. Prior to inoculation they were streaked 
nucleotides from the gene of interest [51]. Both reverse on Luria Bertani agar plates and incubated for one day at 


Table 6 Levels of the NAD + , NADH, NADP + , and NADPH cofactors and their ratios of various P. putida strains under 
PHA producing conditions 


Strain 




Cofactor level (mmol/g rCDW)* 


Cofactor ratio 




NAD + 


NADH 


NADP + NADPH 


NADH/NAD + 


NADPH/NADP + 


KT2440 


7.62 ± 0.86 


0.67 ± 0.08 


3.22 ± 0.38 2.67 ± 0.07 


0.09 


0.83 


KX2440Agcc/ 


6.01 ± 0.67 


1.75 ±0.36 


4.49 ±0.1 2 1.31 ±0.32 


0.30 


0.29 


\<J2440-acoA 


6.25 ± 0.42 


2.02 ±0.1 5 


6.48 ± 0.48 0.87 ± 0.27 


0.32 


0.13 


KY2440Agcd-ocoA 


6.88 ± 0.33 


2.94 ±0.1 5 


3.98 ±0.46 0.96 ±0.36 


0.43 


0.24 


KT2440-pg/ 


5.59 ±0.14 


1.82 ±0.10 


2.58 ± 0.40 3.07 ± 0.41 


0.33 


1.19 


KJ2440 Agcd-pgl 


10.81 ±0.50 


2.1 8 ±0.38 


2.30 ±0.1 5 3.71 ±1.23 


0.20 


1.62 



*rCDW: residual cell dry weight. 
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30°C. Single colonies were then picked from the plate and 
inoculated in 50 mL shake flasks containing 10 ml LB and 
incubated overnight under aerobic conditions at 30°C and 
shaking at 180 rpm. To start the PHA-producing process 
1 L baffled shake flasks containing 200 mL of defined 
culture medium were inoculated to an initial OD 595 of 
0.05 and placed in a rotary shaker (180 rpm) under 
aerobic conditions at 30°C Each culture was carried out 
by triplicate. P. putida strains were grown in a defined 
M9 mineral medium consisting of per liter 12.8 g 
Na 2 HP0 4 7H 2 0, 3 g KH 2 0 4 , 1 g NH 4 C1, and 0.5 g 
NaCl. This basic solution was autoclaved and subsequently 
supplemented with 0.12 g of MgS0 4 H 2 0 and a trace 
element solution made of: 6.0 mg/L FeS0 4 7H 2 0, 2.7 mg/L 
CaC0 3 , 2.0 mg/L ZnS0 4 H 2 0, 1.16 mg/L MnS0 4 H 2 0, 
0.37 mg/L CoS0 4 7H 2 0, 0.33 mg/L CuS0 4 5H 2 0, 0.08 mg/L 
H3BO3, and 20 g/L glucose (all filter-sterilized). Antibiotics 
were 62.5 [ig/mL streptomycin (Sm) for E. coli XLIO-Gold, 
100 (ig/mL streptomycin for P. putida strains and 
100 ug/mL ampicillin (Ap) for R coli DHlOb. 

Bioreactor fermentations 

P. putida strains were grown in M9 medium supplemented 
with 18.5 g/L glucose, and antibiotic as required. Bioreactor 
batch fermentations were carried out in a 2 L top-bench 
BIOSTAT Bl bioreactor (Sartorius B Systems GmbH, 
Melsungen, Germany) with a working volume of 1.5 L, at 
30°C. The aeration rate was set to 500 mL L" 1 min" 1 using 
a mass flow controller (PR4000, MKS Instruments, 
Wilmington, MA, USA). The dissolved oxygen level 
was kept above 20% air saturation by control of the 
agitation speed up to a maximum of 900 rpm. The 
pH was maintained at 7.0 by automatic pH controlled 
addition of 500 mM H 2 S0 4 or 1 M of KOH. 

Cell growth was monitored via OD measurements at 
600 nm (Ultraspec 2000, Hitachi, Japan). The cell dry 
weight was determined gravimetrically after harvesting 
bacteria from 10 mL culture broth for 10 min at 4°C and 
10,000 rpm (Eppendorf 5810 R, Hamburg, Germany), 
washing with distilled water, and drying of the obtained 
pellet at 100°C. The ammonium concentration in cell-free 
extracts was measured by a photometric test (LCK 303 
kit, Hach Lange, Danaher, USA). The glucose concentra- 
tion in the cultivation supernatant was analyzed after 
appropriate dilution by HPLC chromatography (HPLC 
Agilent 1260, Agilent, Krefeld, Germany) through an 
8 mm Rezex ROA-organic acid H column (Phenomenex, 
USA) at a flow rate of 0.5 mL-min" 1 at 65°C using 
0.013 N H 2 S0 4 as the mobile phase, followed by detection 
using a RID detector (Agilent series 1260). Gluconate and 
2-keto-gluconate were quantified by HPLC chromatog- 
raphy using an Aminex HPX 87 H column (Biorad, 
Hercules, CA, USA) on a Hitachi HPLC system (Hitachi 
Elite LaChrom, Krefeld, Germany) and 12.5 mM H 2 S0 4 as 



the mobile phase at a flow rate of 0.5 mL min 1 at 45°C. 
Compounds were quantified by UV detection at 220 nm. 

PHA characterization and quantification 

Composition analyses of the polymer produced and cellular 
PHA content concentration determination were performed 
by gas chromatography (GC) and mass spectrometry (MS) 
of the methanolyzed polyester. For PHA production corre- 
sponding P. putida strains were streaked on LB agar plates 
and incubated for one day at 30°C. Single colonies were 
then picked from the plate, inoculated in a 50 mL shake 
flask containing 10 mL of the above described medium and 
incubated overnight under aerobic conditions at 30°C and 
180 rpm (Innova, Enfield, USA). To begin the PHA- 
accumulating process, 1 L baffled shake flasks containing 
200 mL of culture medium were inoculated at an initial 
OD 595nm of 0.05 and incubated under aerobic conditions at 
30°C. Subsequently, 10 mL of the culture were centrifuged 
for 10 min at 4°C and 9,000 x g, followed by a washing step 
with distilled water and lyophilization. Methanolysis was 
carried out by suspending 5-10 mg of lyophilized cell pel- 
let in 2 mL of chloroform and 2 mL of methanol containing 
15% sulfuric acid and 0.5 mg • mL" 1 3-methylbenzoic acid 
as internal standard. Incubation followed at 100°C for 4 h. 
After cooling, 1 mL of demineralized water was added. The 
organic phase containing the resulting methyl esters of the 
monomers of the polymers of interest were analyzed by 
GC-MS on a Varian GC-MS system 450GC/240MS ion 
trap mass spectrometer operated by the software MS 
Workstation 6.9.3 (Varian Inc., Agilent Technologies, 
Boblingen, Germany). For this purpose, 1 mL of the organic 
phase was injected into the gas chromatograph at a split ra- 
tio of 1:10. Separation of the analytes of interest (i.e. the 
methyl esters of 3-hydroxyexanoate, 3-hydroxyoctanoate, 
3-hydroxydecanoate, 3-hydroxydodecanoate, 3-hydroxy-5- 
cis-dodecanoate, 3-hydroxytetradecanoate) was achieved 
by chromatography through a FactorFour VF-5 ms capil- 
lary column (30 m X 0.25 mm i.d. X 0.25 mm film thick- 
ness). Helium was used as carrier gas at a flow rate of 
0.9 mL • min" 1 . The injector and transfer line temperature 
were 275°C and 300°C, respectively. The temperature pro- 
gram was: initial temperature 40°C for 2 min, then from 
40°C up to 150°C at a rate of 5°C min" 1 and finally up 
to 280°C at a rate of 10°C min" 1 . Positive ions were 
obtained using electron ionization at 70 eV. The mass spec- 
tra were generated by scanning ions of m/z 50 to m/z 650. 
The PHA content (%wt) was defined as the percentage of 
polyhydroxyalkanoate in the cell dry weight (CDW). 

NADH and NADPH quantification 

The intracellular NAD + , NADH, NADP + and NADPH con- 
centrations were measured employing an NADP/NADPH 
Colorimetric Quantitation Kit (Biovision, Abeam, UK). 
Crude cell lysates were obtained by undertaking various 
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freeze-thaw lysis cycles. Cofactors concentrations were 
photometrically quantified by measuring colorimetric 
changes at 450 nm. 

Transcriptome analysis via RNA seq 

Isolation of total RNA from corresponding 10 mL P. 
putida culture was performed using the RNeasy kit 
(Qiagen, Venlo, Netherland), according to the instructions 
provided by the manufacturer. Library preparation was 
done as described in [51]. Briefly, ribosomal RNA removal 
was performed with the MICROExpress Bacterial RNA 
enrichment kit using the Pseudomonas module (Ambion, 
Life Technologies, Darmstadt, Germany) according to the 
manufacturers instructions. The RNA was fragmented to 
a median fragment size of 200 nt by sonication using 
Covaris Adaptive Focused Acoustics device (Covaris, LGC 
Genomics, Herts, UKStrand-specific cDNA libraries 
were generated by previous RNA-Adapter ligation and 
subsequent reverse transcription. Hereby the ligated 
5' -RNA- Adapter contained one of 24 different 6 nt 
barcode sequences [51], which allowed multiplexing of 
several libraries in a single Alumina lane for sequencing. 
Up to 20 libraries were pooled and treated with 
duplex-specific nuclease (DSN) for additional rRNA 
removal. Sequencing was performed on an Illumina 
HiSeq generating single-end reads with a length of 50 
base pairs. Computational analysis of the raw sequence 
data obtained in Illumina FASTQ-format was performed 
as described previously by [53], with some modifications 
Mapping was performed using stampy, a short-read aligner 
that allows for gapped alignments [54] for quantification of 
gene expression. The reads per gene (rpg) values of all 
genes were calculated from the SAM output files. Testing 
for differential expression against KT2440 (four biological 
replicates) was performed with DESeq [55], an R software 
package that uses a statistical model based on the negative 
binomial distribution. 

Protein quantification 

After 60 hours, 10 mL cultures were withdrawn and cells 
were harvested by centrifugation. Proteins from up to 
50 \ig cell pellet were separated via SDS-PAGE [56]. 
Resulting SDS gels were stained using Coomassie Brillant 
Blue and protein amounts quantified densiometrically on 
a GS-800 scanner device. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

IPC designed the study and composed the manuscript. JB and IPC have 
made substantial contributions to the acquisition and analysis of data. AB 
and SH have made contribution of acquisition of the transcriptome data. MS 
and CW have been involved in the analysis of the data. MJ, JB, IPC, and DJ 
have been involved in writing the manuscript critically for important 
intellectual content. All authors read and approved the final version of the 
manuscript. 



Acknowledgements 

We would like to thank V.de Lorenzo and K. Timmis for the constant 
intellectual advice. We acknowledge Esteban Martinez for providing the 
plasmids used in this study. We also thank Monika Schniederjans for the 
excellent support during the transcriptome analysis. We gratefully thank 
SFB-TR51 of the Deutsche Forschungsgemeinschaft for funding the project. 
Ignacio Poblete-Castro acknowledge financial support by the program 
"Convenio de Desempeho Apoyo a la Innovacion en Educacion Superior" 
(PMIUAB 1301). Christoph Wittmann acknowledges financial support by the 
German Federal Ministry of Education (BMBF) via the project "Bio-Nylon" 
(No. 03 V05757). 

Author details 

institute of Microbiology, Technische Universitat Braunschweig, D-38106 
Braunschweig, Germany, department of Molecular Bacteriology, Helmholtz 
Centre for Infection Research, D-38124 Braunschweig, Germany, institute of 
Systems Biotechnology, Saarland University, D-66123 Saarbrucken, Germany. 
4 Universidad Andres Bello, Facultad de Ciencias Biologicas, Biosystems 
Engineering group, Santiago 8340176, Chile. 5 Microbial Drugs group, 
Helmholtz Centre for Infection Research, D-38124 Braunschweig, Germany. 

Received: 17 April 2014 Accepted: 6 June 2014 
Published: 19 June 2014 

References 

1. Mukhopadhyay A, Redding AM, Rutherford BJ, Keasling JD: Importance of 
systems biology in engineering microbes for biofuel production. 

Curr Opin Biotechnol 2008, 19(3):228-234. 

2. Lee JW, Na D, Park JM, Lee J, Choi S, Lee SY: Systems metabolic 
engineering of microorganisms for natural and non-natural chemicals. 
Nat Chem Biol 2012, 8(6):536-546. 

3. Otero JM, Nielsen J: Industrial systems biology. Biotechnol Bioeng 2009, 
105(3):439-460. 

4. Stephanopoulos G: Metabolic Fluxes and Metabolic Engineering. 
MetabEng 1999, 1(1):1-11. 

5. Hln C, Mendes F, Vasconcelos I, Francois JM: A metabolic and genomic 
study of engineered Soccharomyces cerevisioe strains for high glycerol 
production. Metab Eng 2007, 9(4):364-378. 

6. Gao X, Chen J-C, Wu Q, Chen G-Q: Polyhydroxyalkanoates as a source of 
chemicals, polymers, and biofuels. Curr Opin Biotechnol 201 1, 22(6):768-774. 

7. Poblete-Castro I, Becker J, Dohnt K, Santos VM, Wittmann C: Industrial 
biotechnology of Pseudomonas putida and related species. Appl Microbiol 
Biotechnol 2012, 93(6):2279-2290. 

8. Madison LL, Huisman GW: Metabolic Engineering of Poly(3-Hydroxyalkanoates): 
From DNA to Plastic. Microbiol Mol Biol Rev 1999, 63(1 ):2 1-53. 

9. Chung A, Liu Q, Ouyang S-P, Wu Q, Chen G-Q: Microbial production of 
3-hydroxydodecanoic acid by pha operon and fadBA knockout mutant 
of Pseudomonas putida KT2442 harboring tesB gene. Appl Microbiol 
Biotechnol 2009, 83(3):5 1 3-5 1 9. 

10. Escapa IF, Morales V, Martino VP, Pollet E, Averous L, Garcia JL, Prieto MA: 
Disruption of (3-oxidation pathway in Pseudomonas putida KT2442 to 
produce new functionalized PHAs with thioester groups. Appl Microbiol 
Biotechnol 201 1, 89(5):1 583-1 598. 

11. Sun Z, Ramsay J, Guay M, Ramsay B: Carbon-limited fed-batch production 
of medium-chain-length polyhydroxyalkanoates from nonanoic acid by 
Pseudomonas putida KT2440. Appl Microbiol Biotechnol 2007, 74(1):69-77. 

12. Hartmann R, Hany R, Pletscher E, Ritter A, Witholt B, Zinn M: Tailor-made 
olefinic medium-chain-length poly[(R)-3-hydroxyalkanoates] by 
Pseudomonas putida GPo1: Batch versus chemostat production. 
Biotechnol Bioeng 2006, 93(4):737-746. 

13. Poblete-Castro I, Escapa I, Jager C, Puchalka J, Chi Lam C, Schomburg D, 
Prieto M, Martins dos Santos V: The metabolic response of P. putida 
KT2442 producing high levels of polyhydroxyalkanoate under 
single- and multiple-nutrient-limited growth: Highlights from a 
multi-level omics approach. Microb Cell Fact 201 2, 1 1 (1):34. 

14. Tripathi L, Wu L-P, Dechuan M, Chen J, Wu Q, Chen G-Q: Pseudomonas 
putida KT2442 as a platform for the biosynthesis of polyhydroxyalkanoates 
with adjustable monomer contents and compositions. Bioresour Technol 201 3, 
142:225-231. 

1 5. Ouyang S-P, Luo RC, Chen S-S, Liu Q, Chung A, Wu Q, Chen G-Q: Production 
of Polyhydroxyalkanoates with High 3-Hydroxydodecanoate Monomer 



Borrero-de Acuna et al. Microbial Cell Factories 2014, 13:88 
http://www.microbialcellfactories.eom/content/13/1/88 



Page 15 of 15 



Content by fadB and fadA Knockout Mutant of Pseudomonas putida 
KT2442. Biomacromolecules 2007, 8(8)2504-251 1 . 

16. Escapa IF, del Cerro C, Garcia JL, Prieto MA: The role of GIpR repressor in 
Pseudomonas putida KT2440 growth and PHA production from glycerol. 
Environ Microbiol 201 2, 1 5(1 ):93-1 1 0. 

1 7. Klinke S, Dauner M, Scott G, Kessler B, Witholt B: Inactivation of Isocitrate 
Lyase Leads to Increased Production of Medium-Chain-Length 
Poly(3-Hydroxyalkanoates) in Pseudomonas putida. Appl Environ 
Microbiol 2000, 66(3):909-913. 

18. Kenny S, Runic JN, Kaminsky W, Woods T, Babu R, OConnor KE: 
Development of a bioprocess to convert PET derived terephthalic acid 
and biodiesel derived glycerol to medium chain length 
polyhydroxyalkanoate. Appl Microbiol Biotechnol 2012, 95(3):623-633. 

1 9. Muhr A, Rechberger EM, Salerno A, Reiterer A, Malli K, Strohmeier K, Schober S, 
Mittelbach M, Koller M: Novel Description of mcl-PHA Biosynthesis by 
Pseudomonas chlororaphis from Animal-Derived Waste. J Biotechnol 2013, 
165(1):45-51. 

20. Agnew DE, Pfleger BF: Synthetic biology strategies for synthesizing 
polyhydroxyalkanoates from unrelated carbon sources. Chem Eng Sci 2012, 

103:58-67. 

21. Poblete-Castro I, Binger D, Rodrigues A, Becker J: Martins dos Santos VAP, 
Wittmann C: In-silico-driven metabolic engineering of Pseudomonas 
putida for enhanced production of poly-hydroxyalkanoates. 
MetobEng 2013, 15:113-123. 

22. Huijberts GN, Eggink G, de Waard P, Huisman GW, Witholt B: Pseudomonas 
putida KT2442 cultivated on glucose accumulates poly(3-hydroxyalkanoates) 
consisting of saturated and unsaturated monomers. Appl Environ 
Microbiol 1992, 58(2):536-544. 

23. Li Z-J, Cai L, Wu Q, Chen G-Q: Overexpression of NAD kinase in recombinant 
Escherichia coli harboring the phbCAB operon improves poly (3-hydroxybuty rate) 
production. Appl Microbiol Biotechnol 2009, 83(5):939-947. 

24. Yim KS, Lee SY, Chang HN: Synthesis of poly-(3-hydroxybutyrate-co-3- 
hydroxyvalerate) by recombinant Escherichia coli. Biotechnol Bioeng 1996, 
49(5):495-503. 

25. Lee SY, Lee KM, Chan HN, Steinbuchel A: Comparison of recombinant Escherichia 
coli strains for synthesis and accumulation of poly-(3-hydroxybutyric acid) and 
morphological changes. Biotechnol Bioeng 1994, 44(1 1):1 337-1 347. 

26. Lim S-J, Jung Y-M, Shin H-D, Lee Y-H: Amplification of the NADPH-related 
genes zwf and gnd for the oddball biosynthesis of PHB in an E. coli 
transformant harboring a cloned phbCAB operon. J Biosci Bioeng 2002, 
93(6):543-549. 

27. Choi J-C, Shin H-D, Lee Y-H: Modulation of 3-hydroxy valerate molar fraction 
in poly(3-hydroxybutyrate-3-hydroxyvalerate) using Ralstonia eutropha 
transformant co-amplifying phbC and NADPH generation-related zwf 
genes. Enzyme Microb Technol 2003, 32(1 ):1 78-185. 

28. Rojo F: Carbon catabolite repression in Pseudomonas : optimizing 
metabolic versatility and interactions with the environment. FEMS Microbiol 
Rev 2010, 34(5):658-684. 

29. Poblete-Castro I, Rodriguez A, Lam CMC, Kessler W: Improved Production 
of Medium-Chain-Length Polyhydroxyalkanoates in Glucose-Based 
Fed-Batch Cultivations of Metabolically Engineered Pseudomonas putida 
Strains. J Microbiol Biotechnol 2014, 24(1):59-69. 

30. Follonier S, Panke S, Zinn M: A reduction in growth rate of Pseudomonas 
putida KT2442 counteracts productivity advances in medium-chain-length 
polyhydroxyalkanoate production from gluconate. Microb Cell Fact 201 1, 
10(1)25. 

31 . Lee G, Na J: Future of microbial polyesters. Microb Cell Fact 201 3, 1 2(1):54. 

32. Choi J, Lee SY: Factors affecting the economics of polyhydroxyalkanoate 
production by bacterial fermentation. Appl Microbiol Biotechnol 1999, 
51(1):13-21. 

33. Oehmen A, Pinto FV, Silva V, Albuquerque MGE, Reis MAM: The impact of 
pH control on the volumetric productivity of mixed culture PHA 
production from fermented molasses. Eng Life Sci 2013, 1 4(2):1 43-1 52. 

34. del Castillo T, Ramos JL, Rodriguez-Herva JJ, Fuhrer T, Sauer U, Duque E: 
Convergent Peripheral Pathways Catalyze Initial Glucose Catabolism in 
Pseudomonas putida: Genomic and Flux Analysis. J Bacteriol 2007, 

1 89(1 4):5 142-5 152. 

35. Daddaoua A, Krell T, Ramos J-L: Regulation of Glucose Metabolism in 
Pseudomonas: The Phosphorylative Branch and Entner-Doudoroff 
Enzymes are regulated by a Repressor Containing a Sugar Isomerase 
Domain. J Biol Chem 2009, 284(32):21 360-21 368. 



36. del Castillo T, Duque E, Ramos JL: A Set of Activators and Repressors 
Control Peripheral Glucose Pathways in Pseudomonas putida To Yield a 
Common Central Intermediate. J Bacteriol 2008, 190(7):233 1-2339. 

37. Kim J, Jeon CO, Park W: Dual regulation of zwf- 1 by both 2-keto-3-deoxy- 
6-phosphogluconate and oxidative stress in Pseudomonas putida. 
Microbiology 2008, 1 54(1 2)3905-391 6. 

38. Garcia-Mauriho SM, Perez-Martinez I, Amador CI, Canosa I, Santero E: 
Transcriptional activation of the CrcZ and CrcY regulatory RNAs by the CbrB 
response regulator in Pseudomonas putida. Mol Microbiol 2013, 89(1 ):1 89-205. 

39. La Rosa R, de la Peha F, Prieto MA, Rojo F: The Crc protein inhibits the 
production of polyhydroxyalkanoates in Pseudomonas putida under balanced 
carbon/nitrogen growth conditions. Environ Microbiol 2014, 16(1)278-290. 

40. Ryan WJ, O'Leary ND, O'Mahony M, Dobson ADW: GacS-Dependent 
Regulation of Polyhydroxyalkanoate Synthesis in Pseudomonas putida 
CA-3. Appl Environ Microbiol 2013, 79(6):1 795-1 802. 

41 . Berrios-Rivera SJ, Sanchez AM, Bennett GN, San KY: Effect of different levels of 
NADH availability on metabolite distribution in Escherichia coli fermentation 
in minimal and complex media. Appl Microbiol Biotechnol 2004, 65(4):426-432. 

42. Sanchez AM, Bennett GN, San K-Y: Effect of different levels of NADH availability 
on metabolic fluxes of Escherichia coli chemostat cultures in defined medium. 
J Biotechnol 2005, 1 1 7(4):395-405. 

43. Ren Q, de Roo G, Ruth K, Witholt B, Zinn M, Thony-Meyer L: Simultaneous 
Accumulation and Degradation of Polyhydroxyalkanoates: Futile Cycle or 
Clever Regulation? Biomacromolecules 2009, 10(4):9 16-922. 

44. Lee IY, Kim MK, Park YH, Lee SY: Regulatory effects of cellular nicotinamide 
nucleotides and enzyme activities on poly(3-hydroxybutyrate) synthesis in 
recombinant Escherichia coli. Biotechnol Bioeng 1996, 52(6):707-712. 

45. Escapa IF, Garcia JL, Biihler B, Blank LM, Prieto MA: The 
polyhydroxyalkanoate metabolism controls carbon and energy spillage 
in Pseudomonas putida. Environ Microbiol 2012, 14(4):1 049-1 063. 

46. Lee SJ, Song H, Lee SY: Genome-Based Metabolic Engineering of 
Mannheimia succiniciproducens for Succinic Acid Production. Appl Environ 
Microbiol 2006, 72(3):1 939-1 948. 

47. Yim H, Haselbeck R, Niu W, Pujol-Baxley C, Burgard A, Boldt J, Khandurina J, 
Trawick JD, Osterhout RE, Stephen R, Estadilla J,Teisan S, Schrever HB, Andrae S, 
Yang TH, Lee SY, Burk MJ, Van Dien S: Metabolic engineering of Escherichia coli 
for direct production of 1 ,4-butanediol. Nat Chem Biol 201 1, 7(7)445-452. 

48. Bro C, Regenberg B, Forster J, Nielsen J: In silico aided metabolic 
engineering of Saccharomyces cerevisiae for improved bioethanol 
production. Metab Eng 2006, 8(2):102-1 11. 

49. Becker J, Zelder O, Hafner S, Schroder H, Wittmann C: From zero to 

hero - Design-based systems metabolic engineering of Corynebacterium 
glutamicum for l-lysine production. Metab Eng 201 1, 1 3(2):1 59-1 68. 

50. Park JH, Lee KH, Kim TY, Lee SY: Metabolic engineering of Escherichia coli for 
the production of l-valine based on transcriptome analysis and in silico 
gene knockout simulation. Proc Natl Acad Sci 2007, 1 04(1 9):7797-7802. 

51. Silva-Rocha R, Martinez-Garcia E, Calles B, Chavarria M, Arce-Rodnguez A, 
De las Heras A, Paez-Espino AD, Durante-Rodnguez G, Kim J, Nikel PI, Platero R, 
De Lorenzo V: The Standard European Vector Architecture (SEVA): a 
coherent platform for the analysis and deployment of complex prokaryotic 
phenotypes. Nucleic Acids Res 201 2, 41 (D1 ):D666-D675. 

52. Huff JP, Grant BJ, Penning CA, Sullivan KF: Optimization of routine 
transformation of Escherichia coli with plasmid DNA. BioTechniques 1990, 
9(5)570-572. 574, 576-577. 

53. Dotsch A, Eckweiler D, Schniederjans M, Zimmermann A, Jensen V, Scharfe M, 
Geffers R, Haussler S: The Pseudomonas aeruginosa Transcriptome in 
Planktonic Cultures and Static Biofilms Using RNA Sequencing. PLoS One 2012, 
7(2):e31092. 

54. Lunter G, Goodson M: Stampy: a statistical algorithm for sensitive and fast 
mapping of lllumina sequence reads. Genome Res 201 1, 21 (6):936-939. 

55. Anders S, Huber W: Differential expression analysis for sequence count 
data. Genome Biol 2010, 11(10):R106. 

56. Righetti PG: Recent developments in electrophoretic methods. J Chromatogr 
A 1990,516(1)3-22. 



doi:1 0.1 1 86/1 475-2859-1 3-88 

Cite this article as: Borrero-de Acuha et al.: Production of medium chain 
length polyhydroxyalkanoate in metabolic flux optimized Pseudomonas 
putida. Microbial Cell Factories 2014 13:88. 



